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1. Introduction 

The notion of surface tension, or interfacial free energy per unit area, plays a key 
role in many studies concerning the surface phenomena and the phase coexistence. 

When we consider a solid or a fluid, which is a mixture of two chemical species 
1 and 2, in equilibrium with its vapor, one of the problems, experimentally as well 
as theoretically, is to determine how the corresponding surface tension depends on 
the composition of the mixture. 

Some relationship is expected which would give this surface tension, here denoted 
r (i,2)|o, as an interpolation between the two surface tensions, t\\q and T2\o, of each 
of the species when they are chemically pure. 

Using thermodynamical considerations several equations have been derived in 
the literature, according to different assumptions. 

Thus, for ideal or nearly ideal solutions, a fairly simple treatment, due to Guggen- 
heim |TJ, leads to the following equation 

e -/3aV 2) | = Cl e-^ a2T no + c 2 e-^ (1.1) 

where c\ is the fixed molar fraction of species 1 in the (1,2) mixture, C2 = 1 — ci, 
the fixed molar fraction of species 2, a 2 is the mean surface area per molecule, and 
ft = 1/kT is the inverse temperature. 

A very simple relationship for the so called regular solutions comes from Prigogine 
and Defay [2], who proposed the equation 

T(i,2)|o = CiTi|o + c 2 r 2 |o - Kdc 2 (1.2) 

with K a semiempirical constant. 

A simple treatment due to Eberhart |2j assumes that the surface tension of a 
binary solution is linear in the surface composition, that is 

T(i,2)|o = cln\o + c|r 2 |o (1-3) 

where the c|, i = 1, 2, denote the mole fraction near the surface of phase separation, 
and that the ratio c*/ci is proportional to the ratio c s 2 /c2. 

Finally, when the surface tensions t\\q and r 2 |o differ appreciably, a semiempirical 
equation attributed to Szyszkowsky (|3], [H]) gives: 

T(1 ' 2)|0 -l-filnfl + ^ (1.4) 



Ti|o ^ A 

where two characteristic constants A and B of the compounds have been used, and 
c 2 is the concentration of the species with the smaller surface tension. 

We refer the reader to Adamson's book |B| (Chapter III, Section 4), and references 
therein, for a detailed discussion of the above equations. On the other hand, an 
extensive development for various types of non ideal solutions that has been made 
by Defay, Prigogine and co-workers, can be found in the monography 0. 



ON THE STATISTICAL MECHANICS AND SURFACE TENSIONS OF BINARY MIXTURES 3 



More recently, an interface model with a two-valued random interaction was 
introduced by two of the present authors in ref. |B| to describe the phase boundary 
from a microscopic point of view. The surface tension for that model could be 
computed according to a quenched or annealed disorder and one obtains 

quenched 
T (l,2),0 = C 1 T 1,0 + C 2T2,0 

/3„_ annealed Q n 

e -^(i, 3 )|o = cie~ 0T1 '° + c 2 e- l3T2 -° 

in agreement with the above equations or (|1.3j) . 

The aim of the present paper is two discuss the problem within a lattice bulk 
statistical mechanical model describing the binary mixture in equilibrium with its 
vapor. Previous studies of various models of binary lattice gases can be found in 
Refs. Eg. 

Here, we consider a lattice gas system with two kinds of particles, where each 
lattice site can be in one of the three states, 0, 1,2, interpreted, respectively, as an 
empty site, a site occupied by a particle of the first kind of the model, and a site 
occupied by a particle of the second kind. Whenever the particles 2 are not allowed 
the system reduces to the usual Ising model, in its lattice gas version, with coupling 
constant J\/2. We consider the system in the phase coexistence region and denote 
by Ti[o the corresponding surface tension between the dense and the dilute phases. 
Analogously, when particles 1 are not allowed, it reduces to the Ising model with 
coupling constant J2/2 and we let r 2 | be the corresponding surface tension. 

We can also study our three state model in the phase coexistence region (with 
the help of Pirogov Sinai theory) and then interpret the dense phase as the binary 
mixture, the dilute phase as the corresponding vapor, and T(i,2)|o as the surface 
tension between these two phases. On the other hand the concentration of particles 
1 and 2 in the dense phase can be fixed to take any given values. 

As a main result of this paper we prove that, at low temperatures, the following 
equation holds, for the surface tension of our model, 

e -/3(r ( i, 2) |o-^) = c * e -/3(-i|0-^i) + c * e -/% 2 |0-^2) Q 5) 

Here J 7 ^ (i = 1,2) is the specific free energy of the gas of "jumps" describing the 
Gallavotti's line of phase separation for the Ising model in two dimensions [TT] . and 
that of the gas of the "walls" describing the Dobrushin's microscopic interface [12] 
in three dimensions. This means that Ti\ — T\ = J\ and r 2 |o — T2 = Ji are the 
respective energy costs per unit length or unit area of the 1|0 and the 2|0 interfaces. 
The quantity T is the specific free energy (which can be expressed as a convergent 
series at low temperatures) of a gas of some geometrical objects called aggregates. 
In dimension d = 2, those aggregates are the natural generalizations to our model 
of the jumps of Gallavotti's line and the leading term of the series giving this free 
energy T is 

2 c^e- 2/3Jl + cle~ 2Si - h 
~P c\e~^ + c* 2 e-P J z 
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In dimension d — 3, they are the natural generalizations of the walls of the Do- 
brushin's interface and then the leading term of the series is 

1 c\e-^ Jl + c* 2 e~ bl3j2 1 (c^e" 2/3Jl + c^e" 2/3j2 ) 4 
~p c\e~^ + (%e-P J * ~ (djpM + c*e"^) 4 

The coefficients c\ and (f 2 are related to the concentrations c\ and c 2 of the particles 
1 and the particles 2 through equation (|5.7J) below. This equation gives at low 
temperatures: 



(cif 



2dc l e-^(c 1 e-^+c 2 e-^) 2d ~ 1 



2{d + l)ci( Cl e- pJl + c 2 e- /3j2 ) M + 0(e 



-(2d+l)/3min{Ji,J 2 } N 



(1.6) 



for i = 1,2. 



The paper is organized as follows. The model is defined in Section 2 which also 
provides the analysis of the ground states of the system. Section 3 is devoted to the 
study of the Gibbs states of the system at low temperatures and of the coexistence 
between the mixture and the vapor. Section 4 contains the definitions of the surface 
tensions and an expansion of the surface tension between the mixture and the 
vapor in terms of interfaces. Section 5 contains the presentation of the relationship 
between surface tensions. The proofs are given in the two remaining sections. 



2. The model 

We consider a cubic lattice Z d , of dimension d — 2, 3, and to each site Z d we 
associate a variable s x which taking values in the set Q = {0, 1,2} specifies one of 
the possible three states of the system at each lattice site. We say that the site x 
is empty if s x = and that is occupied by a particle of kind 1 or of kind 2 if s x = 1 
or 2. The following Hamiltonian 

2 2 

H = - E E E -A^ oi)5{s y) (3) (2.1) 

(x,y) a=0 /3=0 

where (x, y) denote nearest neighbor pairs, 5 is the usual Kronecker symbol, S(s, s') = 
1 if s — s' and 6(s, s') = otherwise, and E a p = Ep a are the coupling constants, 
is of the form of the Blume-Emery-Griffiths model [12] which describes a general 
three state lattice system for the case of nearest neighbour interactions. We shall 
assume here that 

2E\2 = En + -^22 (2-2) 

in order to ensure that particles of kinds 1 and 2 could be mixed arbitrarily without 
any cost of energy. When these identities are not satisfied two new thermodynamic 
phases, either rich in particles of kind 1 or in particles of kind 2, may appear as 
equilibrium states of the system. 
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With the assumption of hypothesis 1)2.21) the general Hamiltonian (12.11) can be 
reduced to the form 



{x,y) 



Ji(8(s x ,l)5(s y ,0) + 8(s x ,0)5(s y ,l)) 



+J 2 (S(s x , 2)5{s y , 0) + 5{s x , 0)S{s v , 2)) 



(2.3) 



that is, the case in which the coupling constants satisfy E w = E m — —J\. E 20 = 
E 02 = — J 2 and E a p = otherwise. Furthermore, we assume that J\ and J 2 are 
positive constants. 

In order to see this fact we consider, as it is often convenient, the reformulation 
of the three state lattice system in the language of a magnetic system of spin one. 
To do so, we define the spin variable a x at the x site via 

6(s x ,0) = l-a 2 x 

6(s x ,l) = <r x (a x + l)/2 (2.4) 
5(s x ,2) = a x (a x -l)/2 

so that a x = 0, 1, —1 corresponds to the presence at site x of the state 0, 1 or 2, 
respectively. In terms of the spins, the general Hamiltonian 1)2.1)1 takes the form 



-J2(Aa x + Bal) (2.5) 

X 

The pair interactions of the system are given by 

En + E 22 — 2Ei 2 = & J 

En + E 00 - 2E i = 2 J + K + 2C (2.6) 
E 22 + Eoq — 2Eq 2 = 2 J + /C — 2C 

The last two terms can be treated as chemical potentials, with 

A = d(E m - E 02 ) and B = d{E m + E 02 - 2E 00 - 2 J). 

We see that our hypothesis 1)2.2)1 implies 

J = (2.7) 
On the other hand, taking into account that 

25(a x , l)S(a y , 0) = a x (a x + 1)(1 - a 2 y ) = a 2 x + a x - a x a 2 y - o 2 x a 2 
25(a x , -l)6(a y , 0) = a x (a x - 1)(1 - a 2 ) = o 2 x - o x + a x a 2 - a 2 x a 2 

we obtain Hamiltonian 1)2.3)1 . plus chemical potential terms, with 

2Ji = K + 2C, 2 J 2 = K - 2C (2.8) 
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We notice that condition (|2.7jl excludes the models discussed in the above men- 
tioned refs. [10J. In ref. (XQJ , Lebowitz and Gallavotti have considered the cases 
2J = -K > 0, C = 0, also studied by Wheeler and Widom |S|, and J > 0, 
K — C — 0, usually known as the Blume-Capel model. They prove then the ap- 
pearance at low temperatures of two phases, respectively reach in particles of kind 
1 or of kind 2. When the two phases coexist one has, as a consequence of the 
condition C = 0, the equality between the two surface tensions Ti| and r 2 \ . These 
models, therefore, would not be appropriate for the present study. The case C ^ 
(and J > 0) will be briefly commented in Section 5. 

Let us now return to the discussion of Hamiltonian (|2.3j) . Fixed densities of the 
three species are introduced through the canonical Gibbs ensemble of configurations 
sa = {s x } x&A in a finite box A C Z d , such that 

J2 S ( S *>°) = iV o, ^5(s x ,l) = K and ^5(s x ,2) = N 2 (2.9) 

Here No, Ni and N 2 are nonnegative integers satisfying N +Ni+N 2 = |A| where |A| 
denotes the number of sites of A. The associated partition functions with boundary 
condition be are given by 

Z bc (A;N u N 2 ) = J2 e^ H ^6 (^5(^,1),^) 6 ( ^ 5( Sa; , 2), iV 2 ) X bc (s A ) 

(2-10) 

where H\{s A ) is the Hamiltonian (|2.3jl with the sum over nearest neighbours pair 
d(x, y) C A and x bc ( s A) is a characteristic function standing for the boundary con- 
dition be. We shall be interested in particular to the following boundary conditions: 

• the empty boundary condition: X emp ( s A) = FLedA $(. s xi 0) 

• the mixture boundary condition: x mixt ( s A) = ELedA^ ~ ^( s x,0)) 

• the free boundary condition: X fr ( s A) = 1 

Hereafter, the boundary dA of the box A is the set of sites of A that have a nearest 
neighbour in A c = 7L d \ A. 

We define the free energy per site corresponding to the above ensemble as a 
function of the densities p\ and p 2 of the particles 1 and 2: 

f(pi iP2 )= lim--i-ln^ c (A;[ Pl |A|],[p 2 |A|]) (2.11) 
AT^ d p|A| 

where [ • ] denotes the integer part and the thermodynamic limit A j 7L d is taken in 
the van Hoove sense [LI] . 

We introduce also a grand canonical Gibbs ensemble, which is conjugate to the 
previous ensemble, and whose partition function, in the box A is given by 

s A e!) A 
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where the real numbers p,\ and p 2 replace as thermodynamic parameters the den- 
sities pi and p 2 . We define the corresponding specific free energy, the pressure, as 
the limit 

p(p 1 ,p 2 ) = lim T-rrln H 6c (A;//i,/i 2 ) (2.13) 
ATZ d | TV | 

The equivalence of the two above ensembles is expressed in the following 

Theorem 1. Limits ^2.11\) and ^2.1^1 . which define the above free energies, ex- 
ist. They are convex functions of their parameters and are related by the Legendre 
transformations 

p(pi,p 2 ) = sup[/iipi + p 2 p2 - I3f(pi,p 2 )] (2.14) 
Pf(pi, P2) = sup [pxpx + P2P2 - p(pi,p 2 )] (2.15) 

Proof. We consider two parallelepipedic boxes A' and A" of both the same size and 
paste them to form a parallelepipedic box A = A'UA" in such a way that A' PI A" = 
and each site of some side of A' is a nearest neighbour of a site of a side of A". It 
is easy to see that the following sub-additivity property holds: 

Z emp (A; N[ + N», N' 2 + N») > Z emp (A'; N[, A^)Z emp (A"; N», N») 

The same property is shared by the partition function with mixt boundary condi- 
tions. Then the statements of the theorem follow from standard arguments in the 
theory of the thermodynamic limit ^1]. □ 

We next introduce the finite volume Gibbs measures (a specification) associated 
with the second ensemble: 

p -PH A (s A ) bc(„ \ 

P A c (s A ) = 6 „ X l ' Aj (2.16) 
^ bc (A; pi, p 2 ) 

where 

^ A (s A ) = F A (s A ) -jYl 6 ^ V-jT, 6 ^ 2 ) ( 2 - 17 ) 

They determine by the Dobrushin-Landford-Ruelle equations the set of Gibbs 
states Gp(H) on Z d corresponding to the Hamiltonian H at inverse temperature (5 
(see e.g. |15f). If a Gibbs state P G Qp{H) happens to equal the limit lim A | Z d P A C (-), 
we shall call it the Gibbs state with boundary condition be. 

In the zero temperature limit the Gibbs state with empty boundary condition is 
concentrated on the configuration with empty sites: 

lim P A c (emp A ) = 1 (2.18) 

p— >oc 

where emp A is the configuration where all the sites of A are empty, and this limit 
vanishes for any other configuration. Gibbs states at (3 = 00 will be called ground 
states. 
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Let 

Rf^ = [ s g Q A : Vx G A, s x + 0} (2.19) 

be the restricted ensemble of configurations in A with non empty sites, and _R A lixt (c), 
< c < 1 the subset of configurations of _R A lixt with exactly [c|A|] = N sites occupied 
by a particle of the specie 1 (and |A| — [c|A|] sites occupied by a particle of the specie 
2). The number of configurations i? A ixt (c) equals the binomial coefficient ('^')- 

With the mixture boundary conditions one has 

Mi[c|A|] M 2[(l-c)|A|] 

lim Pf xt (s A ) = —. nrr- for each s A e i?f xt (c) (2.20) 

/3^oo A v ' ( e Mi+e^ 2 )l A l v ' 

and 

(W\ el MiN e H2(\A\-N) 

lim Ff^HRrHc)) = ^n— (2-21) 

while this limit vanishes for those Sa G" i? A ixt (the denominator in (I2.20JI and (I2.21|) 
is the sum ^^Lq of the numerator of the R.H.S. of (|2.21jl ). Notice that all config- 
urations sa G i? A lxt (c) have the same probability. Moreover, by Stirling's approxi- 

|A| 

, and the maximum of 



mation one has for large |A| that (JJm) ~ {^Y (l^) 
e^ ic /(e m + e M2 ) is reached for 



1-c 



(2.22) 



The principle of maximal term gives that, for such values, (I2.21jl tends to 1 in 
the thermodynamic limit. This means that the ground state with mixt boundary 
conditions is concentrated on the restricted ensemble i? mixt (c) of configurations of 
non empty sites with concentration c of particles 1 and concentration 1 — c of 
particles 2. 

With free boundary conditions, one has 

lim P A r (emp A ) = -. (2.23) 

AV LAJ i + ( e m+ e ^)l A l v ' 

( |A| )e#*ilA| e /«[(l-c)lA|] 

lim P A r (#? ixt (c)) = ™^ (2.24) 

/3^oo 1 + ( e w + e^ 2 )l A l v 7 

Thus, with the above considerations, we get that for e M1 +e M2 = 1, the configuration 
with empty sites coexists with the restricted ensemble i? mixt (c). The diagram of 
ground states is shown in Figure 1. 
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empty 
configuration 



Figure 1: The diagram of ground states. 
In the next section, we extend this analysis to low temperatures. 



3. Coexistence between the mixture and the vapor 

To extend the analysis of the previous section to the Gibbs states at low tempera- 
tures, we will express the partition functions (|2.12l) with empty and mixt boundary 
condition in term of contour models. 

Let us first introduce the notions of contours by the following definitions. 

Consider a configuration sa with empty sites on the boundary dA (s x = for 
all x G dA). We define the boundary B(sa) as the set of pairs {s x ,s y } such that 
7^ and s y = 0. To a nearest neighbour pair (x, y) let us associate 

a) : in dimension 2, the unit bond (dual bond) b xy that intersects the bond xy 
in its middle and orthogonal to xy. 

b) : in dimension 3, the unit square (dual plaquette) p xy that intersects the 
bond xy in its middle and orthogonal to xy. 

Two pairs {s x , s y } and {s z , s t } of -B(sa) are said adjacent if one of the two conditions 
is fulfilled 

i) : the dual bonds b xy and b zt (respectively the plaquettes p xy and p zt ) are 
connected. 

ii) : x = z, s x = s z ^ 0, s y = s t = 0, and the bond xy with endpoints x and y 
is parallel to the bond with endpoints z and t. 
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A subset B of B{s\) is called connected if the graph that joins all adjacent pairs 
of B is connected. It is called contour of the configuration sa if it is a maximal 
connected component of -B(sa). 

The boundary and the contours of a configuration Sa with occupied sites on the 
boundary of A are defined in the same way. 

A set T of pairs {s x ,s y } such that s x ^ and s y = is called contour if there 
exists a configuration sa such that T is a contour of s\. We use ^(r) to denote 
the set of sites for which s x = a. The set supp T = Sq(T) U S'i(r) U S^r) is called 
support of the contour T. We will also use S(T) = Si(T) U S^fT) to denote the set 
of occupied sites of the contour, Li(T) (respectively L 2 (T)) to denote the number of 
nearest neighbour pairs (x, y) such that such that s x = 1 and s y = (respectively 
s x = 2 and s y = 0) and L{T) = L^T) + L 2 {T). 
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Figure 2: A configuration with two contours: the two rectangles on the left belong to the same 

contour due to the condition ii) of adjacency. 

Consider the configuration sa having T as unique contour. The difference A\5'(r) 
splits in components (set of sites for which the graph that joins all nearest neighbour 
pairs is connected) with either all occupied sites or all empty sites. The component 
that contains <9A, denoted Ext a T is called exterior of the contour. When sa has 
empty sites (respectively occupied sites) on its boundary, T is called emp-contour 
(respectively mixt-contour). The interior of the contour is the set Int T = A \ 
(S'(r) U ExtA T). It is the union of the components of Int emp T with empty sites 
with the components of lnt m i Xt T with occupied sites. Finally V(T) = S(T) U Int T. 
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Two contours Ti and T 2 are said compatible their union is not connected. They 
are mutually compatible external contours if furthermore V(ri) C ExtA T 2 and 
V(T 2 ) C Ext A IY 

With these definitions one gets the following expansions of the grand canonical 
partition functions 

n 

"emp(A; /il, /i 2 ) = |^ ^(rj)"emp(Int e mp Ti! /^l, A t 2)^mixt(I n tmixt Tj; /Xl, /i2) 

{r l5 ...,r„} ext i=i 

(3.1) 

S mixt (A,/z 1 ,/i 2 )= Yl (e^ + en lA ^ V ™ 
{ri,...,r„} ext 

n 

- 

1=1 



J^J ^(rj)r, em p(Int em p Ti] /ii, /i 2 )^mixt(I n tmixt ^i] fJ>l, ^2) (3.2) 



where the first sum is over families of mutually external emp-contours, the second 
sum is over families of mutually external mixt-contours, and 

u(T) = e- /3Jlil ( r )- /3j2i2 ( r ) + ^ 1 l 5l ( r )l + ^ 2 l 5 ' 2 ( r )l (3.3) 

We put 

9emp = 0, fi-mixt = hi(e m + e M2 ), g max = max{g emp , g mixt } (3.4) 

We divide in (|3.1j) each H mixt by S emp and multiply it back again in the form 
(j3.1jh Continuing this process, and doing an equivalent procedure with (13 .2j) . these 
relations lead to the following expansion for the partition functions with boundary 
condition q = emp or q = mixt . 

n 

S 9 (A;/i 1 , / i 2 )=e^l A l U z ^) ( 3 - 5 ) 

{ri,...,r„} corop »=1 

where the sum is now over families of compatible g-contours and the activities z q (T) 
of contours are given by 

z q (F) = ^(De^'^' Tll 11 ^^ 1 '^ (3.6) 

r, q (mtg T; fix, ^2) 

where m 7^ q. 

The (generalized) Peierls estimates 

w ( r)e - 9max |s(r)| < e -/3JL(r) (3 _ 7) 

where J = min{Ji, J 2 }, allow us to have a good control of the behavior of our 
system at low enough temperatures using Pirogov-Sinai theory [IB]. Choosing the 
Zahradnik's formulation of that theory [TJj, we introduce to state our results, the 
following 



p 9 (Hu H2) = lim yyy In E' q (A; pi, p 2 ) (3.10) 
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Definition 1. For q = emp and q = mixt we define the truncated activity 

Z '(T) = { Z ^ if ^( r )<e- Qi(r) 
qV ' \ e ~ aL ^ otherwise 

where a is some positive parameter to be chosen later (see Theorem^) 

Definition 2. The q-contour T is called stable if 

z q (T) < e" aL(r) (3.8) 

i.e. ifz q (T) = z' q (T). 

We define the truncated partition function E' q (A) as the partition function ob- 
tained from (j3.5j) by leaving out unstable contours, namely 

n 

E' q (A; p,, p 2 ) = £ j]^) (3.9) 

f r l r »}comp i = 1 

Here the sum goes over compatible families of stable q -contours. Let 

1 

A(z d \A\ 

be the meta-stable pressure associated with the truncated partition function 
E' g {A;pi,p 2 )- 

For a large enough the thermodynamic limit (13.101) can be controlled by a con- 
vergent cluster expansion. 

Namely, to exponentiate the truncated partition function, we introduce multi- 
indexes X as functions from the set of contours into the set of non negative in- 
tegers (see [18j,[19j). We let supp X = Ur ; x(r)>isupp T and define the truncated 
functional 

n^n%(r)* (r) (3.11) 

where the factor a(X) is a combinatoric factor defined in terms of the connectivity 
properties of the graph G(X) with vertices corresponding to T with X(T) > 1 
(there are X(T) vertices for each such T) that are connected by an edge whenever the 
corresponding contours are incompatible). Namely, a(X) = and hence $ q (X) = 
unless G(X) is a connected graph in which case X is called a cluster and 

a(X) = (-l) |e(G)l (3.12) 

GCG(X) 

Here the sum goes over connected subgraphs G whose vertices coincide with the 
vertices of G(X) and |e(G)| is the number of edges of the graph G. If the cluster 
X contains only one contour, then a(X) = 1. 

Note that the number of contours T with I^COI = s, L(T) = n, whose support 
contains a given site can be bounded by 2 s v§ where z/ 2 = 4 and u 3 = 14?. 
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As a result of standard cluster expansion we get that for ne~ a < 1 where k = 
2i>dK c \ and k c1 = ^ e^ 1 is the cluster constant [2U] : 

lnH;(A;^,// 2 ) = g q \A\ + ^ $ q {X) (3.13) 

X:suppXcA 

= ^iai + iai E j|^ + *( A i**) ( 3 - 14 ) 



X:suppX9z 



where 



E '7 pX "| A k w (3.i5) 

^— ' SUDD A 

X:suppXnAV0 

In addition (see Section 6): 

^ |$ 9 (X)| < Ke~ a (3.16) 

X:suppX9z 

|a(A | $,)| < Kc~ a \dA\ (3.17) 



giving 



Pqim,^) =g g + Yl \^T~x\ ( ^ 3 ' 18 ' ) 

X:supp X3i l SU PP I 

The following theorem shows that the low temperature phase diagram of the 
model is a small perturbation of the diagram of ground states (see Figure 3). 

Theorem 2. Assume [3 is large enough so that e _/3J+5 = e~ a < 2 (d+i)K > then there 
exists a coexistence line ln(e Ml + e M2 ) = 0(e~ 2dl3J ) on which all mixt-contours and 
all emp- contours are stable and such that: 

E 9 (A;^)=E*(A;^) (3.19) 

for both boundary conditions q = mixt and q = emp, and the pressure is given by 

PilH^l) = Pmixt(K, Hi) = Pempi^l, f4) ( 3 - 20 ) 

For any t > 

S mixt (A; fi{ +t,(il + t) = H' mixt (A; \l\ + t,fx* 2 + t), (3.21) 

p(f4 + t, H* 2 + t) = Pmirt(A*l + *, A*2 + *) > Pempi^l + t, f£ + t) 

and 

z>emp(A; /i 1 — t, fi 2 — t) = ^ omp (A, ; fi 1 — t, fi 2 — t), (3.22) 

P{H\ -t,t4-t)= PempilA. -t,t4-t) > Pmixt(^* ~ t, fi* 2 - t) 

The proof is postponed to Section 6. 
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M2A 



mixture 



Ml 



vapor 



Figure 3: Phase diagram at low temperature. 

Whenever a cluster X contains only one contour T (X(T) = 1 and X(T') = for 
I"" 7^ T) one has $ g (X) = z q (T). From this property we get that the metastable 
pressures read: 

p mixt (Mi,M 2 ) = Me* 11 + e^) + 1 +0{e< M+1 » J ) (3.23) 



p emp (Mi,/i 2 ) = e^i-W + e f«-W + o( e -(MWW) 



(3.24) 



Equalizing (|3.23j) with (I3.24JI gives the first term of the equation for the coexistence 
line stated in Theorem [2j 

Let us introduce the infinite volume expectation ( • ) bc (/ii, /i 2 ) associated to the 
Gibbs measure (I2.16J) : 



(■> bc (/xi,/x 2 ) = Jnn Y, -PaCsa) 
As a consequence of the cluster expansion we have for any t > 0: 



(3.25) 



(6(s x ,l) + 5(s x ,2)r^(fi* 1 + t,fi* 2 + t) = (l-5(s x ,0)r^(^+t,fi* 2 + t) 

>l-0(e- 2dl3J ) (3.26) 

(S(s x , 1) + 5(s x , 2)) em P(/il -t,£-t) = (l- 5(s x , 0)) emp (^ - *, /4 - t) 

< 0(e- 2d/3J ) (3.27) 

This shows that the model exhibits at low temperature a first order phase transition 
at the coexistence line where the pressure is discontinuous. 
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Let, Q = ®Es™0l>M) | Mi=M » [At2=M » ; % — i 5 2 be the density of the particle i in the 
mixture regime, on the coexistence line, and let 



-A:supp A9i 



/il=/il,jtl2=/i2 



= 2de^-^ (e ? . + * J 

( e /i* _)_ e^js^+i 



for i — 1, 2. One has 

c < = Wr?f +<i - = ^r^ 1 ^^*' = jra( 1+ °( e " 2 * SJ )) (3 - 28) 

4. Surface tensions 

To introduce the surface tension between the mixture and the vapor, we consider 
the parallelepipedic box: 

V = V LM = {(x!,..,x d ) G Z d : \ Xi \ <L,i = l,...,d-l;-M< x d < M- l} 

Let <9 + V (respectively be the set of sites of 9V with x d > (respectively 

x d < 0). We introduce the boundary condition 

x mixt,e mp(sv0 = Y[ (l-5(s x ,0)) J] 5(s x ,0) 
xed-V xed+v 

This boundary condition enforces the existence of an interface (see below for its 
precise definition) between the mixture and the vapor. The interfacial tension 
between the mixture and the vapor is defined by the limit 

F(V) 

l"^ m"« (2L + 1) 

where 



mixt,e mp = lim lim — — (4.1) 



F(y) (3 ln (s mixt (y ; / ,^ t '^) P s e ^( 1 ^, ^))va (42) 

This definition is justified by noticing that in this expression the volume terms 
proportional to the free energy of the coexisting phases, as well as the terms corre- 
sponding to the boundary effects, cancel and only the term that takes into account 
the free energy of the interface is left. 

To give a precise description of interfaces, we let L + denotes the semi infinite 
lattice with x d > and let L_ = Z d \ L + denotes its complement. Consider then a 
configuration s G ViF 1 , with empty sites on d + V and on L + \ V and with occupied 
sites on a d-V and on L_ \ V. The boundary B(s) (set of pairs {s x , s y } such that 
s x 7^ and s y = 0) of such configuration necessarily contains an infinite component 
loo whose support outside the box V is the set of n.n. pairs between L + and L_. 
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We call interface / the part of Too whose support lies inside the box V: 1^ \ I, is 
called extension of /. As it was done for contours, we use S a (I) to denote the set of 
sites for which s x = a. The set supp / = S (I) U Si(I) U S 2 (I) is called support of 
the interface /. We will also use S(I) = Si(I) U S 2 (I) to denote the set of occupied 
sites of the interface, Li(I) (respectively L 2 {I)) to denote the number of nearest 
neighbour pairs (x, y) such that such that s x = 1 and s y = (respectively s x = 2 
and s y = 0) and L(I) = L^I) + L 2 (I). 
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Figure 4: Interface of a configuration with the mixt, emp boundary condition. 

The set V \ S(I) splits into a part D below the interface and a part U above 
the interface: if one consider the configuration whose boundary contains only the 
interface J, D is the subset of V \ S(I) with occupied sites and U is the subset of 
V \ S(I) with empty sites. 

Then, the partition function, with mixed boundary conditions, can be expanded 
over interfaces as follows 

£mixt,emp(^;/^,/4) = Yl U ( J ) S mbrt (D\ //* , ^)S emp (U] n\ , (4.3) 

/:supp ICV 

where 

W (I) = e -f3JiLi(I)-f3J2L 2 (I)+n* 1 \Si(I)\+ f M* 2 \S 2 (I)\ ( 4 _ 4 ) 

Therefore, 

~/:silcv u (s mixt (y ; ^,^)s emp (y ; ^,^))V2 ^ 
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Since H 9 (A;^*,^|) = S' g (A; /4, /4) f° r both empty and mixt boundary conditions, 
we get by (ETT1 ) and (FTTSl) 

2 p (A; i4, (4) = exp {p q \A\ + a{A \ %)} 

Applying this formula to the various partition functions of (|4.5J) . and taking into 
account equations (|3.20JI . we get: 

e -(3F(V) = ^ CJ (/) e -P(M^^)|S(/)|+a(I)|$ mix 0+a(^|$emp)-| CT (V'|* mixt )-| ^ 7(V'|^.emp) 

-f:supp IC.V 

(4.6) 

Let V$ = {Z d \ V) n L + and Vt = {Z d \ V) n L_. We then apply the formula 
((SUSj) to the four last terms of the exponent of the RHS of (|4.6j) and rewrite the 
different contributions according to appropriated and natural decompositions of the 
sets over which the sums take place and collecting analogous terms. Before applying 
this formula it is convenient first to sum over the clusters with same support. Thus 
we let 

$mixt(C) = Yl $ mixtW 
X:supp X=C 

<W^) = $ emp(X) 
X:supp X=C 

We then get 

e -W) = 6 -k v uj(I)e- p ^^ s ^A(I)B v (I) (4.7) 

/:supp ICV 

where 

A(I)= J] e -*"- t(C)K rSP J] e -*e mp (c)^p (49) 
C:Cns(i)^9 C:Cns(i)^9 

and 

/ x tt j> ■ rr"i I CnD I -n- * . ir\ l Cn ( y \°)l 

S v (7) = ] | e *-xt(C)^- J | e * m , xt (c) — ^ — 

C:CnS(/)^0 C*:CnVf^0 
Cnvf ^0 

gSemplO)^- 1 j e *emp(,Cj p (4.10) 

C:CnS(/)^0 c-.CnV?^® 

In these formula, it is understood that the arguments C of $ m ; x t and $ e m P are 
respectively supports of clusters of mixt-contours and clusters of emp-contours. 
These functionals satisfy the bound: 

|$(C)| < L(C)(Ke- a ) L{c) (4.11) 

where L(C) is the number of n.n. pairs of C (see Section 6). 
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This property allows us to prove that the limit of F(V) when M — ► oo exists, and 
that, if we denote by V the infinite cylinder liniM->oo Vl,m^ ^ nen one S e ts actually 

= F(V) with F(V) defined as in pT^ - fOTljl . 
The surface tension then reads 



F(V) 

W, emp = jhn ( ^ - (4.12) 



Clearly the term Ky/(5{2L + l) d 1 tends to in the limit L — ► oo. Let us 
introduce the modified free energy 

e -/9F'(V) = £ W (/) e -^^)|S(/)U( 7 ) (4 . 13) 
/:supp IC.V 

As a consequence of the analysis of Section 7 we shall see that the free energy 
F'(V) differs from F(V) only by a term proportional to L d ~ 2 , thus showing that 
the surface tension is also given by (I4.12JI with F(V) replaced by F'(V). To simplify 
notations, we shall only consider this last free energy. 

In the Solid-On-Solid (SOS) approximation, that we will also consider, the surface 
tension reads 

T S 9 S = lim — [12 (4 u) 

where 

F sos (F) = ~\n Yl w(/Xe" ! +e"S)- |5(J)l 

i" sos :supp /CV 

where the SOS interfaces belonging to J sos have no overhangs. This means that 
the set of dual bonds (in 2-dimensions) or dual plaquettes (in 3-dimensions) of 
such interface corresponds to the graph of a function. This approximation may be 
obtained by adding to the Hamiltonian (|2.Hjl the anisotropic interaction 

Evert 
J [S(s X} 0)(1 - 6(s v , 0)) + (1 - 6(s x , 0))6(s y , 0)] 

Evert 
is over vertical bonds, and then taking, with an appropriated 

normalization, the limit J 1 — > oo. The coexistence line in that approximation 
coincides with the ground states coexistence line, e^* + = 1, so that: 

e-^ sos ^)= ^( J ) ( 4 - 15 ) 

/ sos :supp ICV 

Let us now define the surface tensions between each species of the mixture and 
the vapor. As mentioned in the introduction, whenever either the particles 1 or the 
particles 2 are not allowed the system reduces to the usual Ising model in its lattice 



ON THE STATISTICAL MECHANICS AND SURFACE TENSIONS OF BINARY MIXTURES 19 




gas version. Thus, we introduce the configurations n v G {0, l} y of the lattice gas 
and the following partition functions 

Qa{V) = ^2 ex P\ P J <* ~ n y) + ( 1 ~ n x 

n v e{o,i} v { {x,y)cv 

Qa,o(V) = ^2 exp<(3J a ^2 [n x (l - n y ) + (I - n x )n y ] 
m,e{o,i} v { (x,y)cv 

x n (i - "•'•» n Ux 

x&d-V xGd+V 

for a = 1 and a = 2. 

The interfacial tension between the species a = 1,2, and the vapor is the limit 

(PI EH) 

where 

F(V)- 1 In Q ^ V) 

It is well known that the ratio Q a fl(V)/Q a (V) can be expressed as a sum over 
interfaces which in this case are connected set of bonds or plaquettes of the dual 
lattice [H1II2]. Extracting the energy of the flat interface, the system can be written 
as a gas of excitations leading to 

F a (V) = J a (2L + l) d - 1 + F^(V) 

In two dimensions, is the free energy of the gas of jumps of the Gallavotti's 
line [TT]. In three dimensions, F^ x is the free energy of the gas of walls of the 
Dobrushin's interface [12] • In both cases these free energies can be analyzed by 
cluster expansion techniques at low temperatures. Namely, the specific free energies 
T a = lim^oo F^ xt (V)/(2L + l)^ -1 exist and are given by convergent expansions in 
term of the activities e~^ Ja , giving 

Pr afi = (3J a + (3T a (4.16) 

In addition 

- $T a = 2e~ f3Ja + 0(e~ 2pJa ) for d = 2 (4.17) 

-PF a = 2e~ ApJa + 0(e- 6/3Ja ) for d = 3 (4.18) 



We refer the reader also to [22], [23j, [24] and references therein for the study 
of these expansions. Furthermore, in two dimensions the surface tension defined 
above is known to coincide with the one computed by Onsager [IB]. We thus have 
an exact expression for r Qj0 , and for T a : 

/3^ a = lntanh(/3J a /2) (4.19) 
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for (3J a larger than the critical value ln(l + y/2). 

Similar results hold in the corresponding SOS approximation (see [25j for the 
three 3-dimensional case). We will use and the surface tensions and 

free energies in this approximation. In two dimensions J-~a° S i s a l so exactly known 
[26] and also given by ()4.19|) . but for all temperatures: 

(3^ os = lntanh(/3J a /2) (4.20) 

for 0>O. 

5. Main Results 

In this section we will give a relationship between the surface tensions introduced 
in the previous section. 

The leading term of the free energy F'(V) corresponds to flat interfaces without 
any decoration. They are those interfaces, for which the set of n.n. x, y such that 
x is empty and y is occupied crosses the plane Xd = —1/2, and such that A(I) = 1. 
Let Xfl at be the set of flat interfaces and N = (2L + l) d_1 . We have 

i:supp iav 

'"flat 



N 'JV 



n 



e (p.i-pj 1 )n e (rt-/3J a )(N-n) e -p(jj,l,i4)N 



where 



n=0 



We will show in Section 7 that the difference F'(V) — F^iV) can be expressed 
as a free energy of a gas of excitations called aggregates. It will then turns out that 
the limit 

T=lim ^n^n (5 . 2) 

exists and is given by a convergent expansion at low temperatures, see (|7.17j) . 

Theorem 3. Assume (3 is large enough so that 8e(e — l)/t 2 e _/3J+5 < 1, then the 
interfacial tensions w^mp, r^o and r 2j o satisfy the equation: 

where 

c* = e ^- p ^*'^ (5.4) 
T is the convergent series defined by \7.17\) : 

-pj: = ° l6 + ° 26 - + 7 6 m + ° 26 R1 . A + 0(e" 5/3J ) (5.5) 
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in dimension d = 3 and 

„* p -2/3Ji , * -2/3J 2 

= 2 Cl " +C f + 0(e" 2 ^) (5.6) 

c * e -/3Ji + c * e -f3J 2 ^ 7 ^ ' 



dimension d = 2, and the convergent series T a satisfies \4.11 ) \4.1£fy and \4.19\) . 



The proof is postponed to Section 7. 
Some remarks and comments are in order. 

The densities q are easily related to the c* through the relation (see 1)3.2 

°* = Q l + (e^ + e^)^ (5 - 7) 



1 _ ( Cie -/9Ji + C2 e-^)^ a - 2rf Q e- M (c 1 e- /3Jl + c 2 e~^) 



(5.J 



In the SOS approximation, as mentioned in the previous section, the coexistence 
line is given by e^* 1 + = 1 and the densities of each species in the mixture regime 
on this line are c\ = e M i and c 2 = e^. We have 

-f:supp ICV 

and the equation between the surface tensions reads 

e -/«JU P -^ sos ) = Cie -/^T-^ sos ) + C2 e-^T-^ os ) (5.9) 

Here jF bUb = liiriL_> 00 N flat : it can be expressed as the series given by 

1)7.1 7)1 . but with truncated functional corresponding to the activities 

z $OS/ w \ _ e -0JiL 1 (W)-/3J 2 L 2 (W)+ t il\S 1 (W)\+^\S2(W)\ 

and satisfy also relations 1)5.5)1 and 1)5.6)1 with c* replaced by Cj. 

To find an exact solution for JF 505 in two dimensions seems to be an interesting 
problem. Actually, this would give an exact equation for the surface tensions not 
restricted to low temperatures [27) . 

The method developed here can be extended naturally to finite range interac- 
tions. As an example, we can consider the Hamiltonian 1)2.3)1 with the sum over 
nearest neighbours and next nearest neighbours. In that case we found that the 
corresponding surface tensions satisfy, in two dimensions and at low enough tem- 
peratures, the equation 

g— ^(Tmixt.emp--^) _ (^g -3 ^ 1 _|_ c\t~^ J2 

— c * e -P(-ri,o-^l) _|_ c * e -P(T2,o-F 2 ) 
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where 



2/3 Ji> 



and 



-I3T 



c f e -W-h + c *2 c */ e -6/3Ji-J 2 + g 



-/3JF 2 = 2e"^ 2 + 0(e- 2 ^ 2 ) 

-5Ji-2J 2 _|_ e -4Ji-3J 2 ^ 



c^ 2 (e" 



3/3Ji-4J 2 _|_ e -2Ji-5J 2 _|_ e -Ji-7J 2 -\ _j_ c *3 e -7/3J 2 



+ 0(e- 



with c* = Q(l + 0(e- 8/3J )). 

When the hypothesis (|2.2jl or (|2.7jl on the model defined by Hamiltonian ()2.1I1 or 
1)2.5)1 is not satisfied, with J > 0, there are, at low temperatures, two phases rich 
in particles of species 1 or of species 2 as been proved in [1QJ for C = 0. This can 
be also shown for C ^ provided that the three quantities in equations 1)2.6)1 are 
positive. Indeed the Hamiltonian 1)2.1)1 and 1)2.5)1 can be reduced, up to chemical 
potential terms, to the form 



H 



E 



J 1 (5(s x ,l)5(s y ,0) + 5(s x ,0)8(s 



y 



1)) 



+ J 2 (5(5*, 2)<f(s tf , 0) + 5(s x , 0)5(s y , 2)) 
+j(<5( Sa ,,l)<y(s v , 2) +5(s x ,2)S(s y ,l)) 



(5.10) 



where Ji = 2J7" + /C + 2C, Ji = 2J + /C — 2C, and J7" are positive constants. In this 
situation and with fixed concentration of the species either we are in one of these 
two pure phases that coexist with the gaseous phase or in the coexistence of the 
three phases. In the later case, there will be segregation between the three phases 
and one will observe only the three kinds of interfaces between these phases. One 
cannot truly speak in this approach about a real surface tension of a mixture. 

Let us stress that from our equations ()5.3)1 and 1)5.9)1 . we get obviously the 
Guggeinheim relation 1)1.1)1 by neglecting the terms exponentially small with the 
inverse temperature, i.e. J 7 ,^, and J-^- This shows that this relation provides a 
good approximation at very low temperatures. 

To see the quantitative difference between the above equations, let us consider a 
concrete example. Note that when we take only into account the (exponentially) 
smallest order corrections in our equations, they reads (in dimension 3): 



T, 



mixt,emp 



kT 

— 5- In 

a 2 



-2e~ 



4 ° 2 n,o 



~TT~ 



+ (l-ci)e" 



T 2,0 



2e~ 



ia 2 T 2 fi 

' — FT — 



(5.11) 



+ kT 



c\e 



+ (1 -ci)e~ 
+ (l-ci)e- 



5 " 2t 2 ,o 



+ kT 



+ (1 -ci)e~ 
+ (l-Ci)e- 



2a T 2,0 
kT 



"I 4 
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where a is the lattice spacing and k the Boltzman constant, while for the Guggen- 
heim relation, one has 



Gugg 



kT 



mixt,emp 



In 



° 2t i,o 



cie ^ + (1 — ci)e 



r 2,0 



(5.12) 



Note also that the simplest form of equations (|1.2jl and t|l . 3|) reads 

T(i,2)|o = cin )0 + (1 - ci)t 2 , (5.13) 

Figure 5 shows that there is actually a quantitative difference between the equations 
1)5.11)) and (|5.12|) . We have chosen a mixture of water and hexane, i.e. r li0 = 

o 

0.0724iV/m, t 2)0 = 0.01978iV/m, a = 3A, the size of a water molecule, and the 
temperature T = 300K. 




Figure 5: Plot of the surface tension given by 15.11|l (solid curve) as function of the concentration 
c\ . The dotted curve corresponds to formula (I5.12|) and the dotted line to formula 15.13|l . 

Finally, notice that when the coupling constants J\ and J 2 (or equivalently the 
surface tensions r^o and r 2i o) differ appreciably, we get assuming J\ < J 2 , and 
neglecting in equations ( 15.11)1 and ( 15.12)1 the terms exponentially small with respect 
to the inverse temperature: 

T mix t,emp =1 _ J_ ln(l _ c *) (514) 
Tl,0 pJl 

to be compared with Szyzkowsky's equation (|1.4)1 . 

6. Proof of Theorem 2 



Let us first give the proof of relations (|3.13)) - (|3.17)1 . 
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Let n(T) = (ave x )- L( - T) , with v = 2v d , a > 1, and A > 0, then 



e A 



£ M(r) < l(t q ) e- Xn *~ n < 73 7 L ( r o) (6.1) 



a 

rvr n=i 



where the sum runs over contours T incompatible with a given contour r . 

The condition ue'^^ae^^ < 1 actually ensures that the convergence condition 



\ Z q 



;(r )| < (e^-l)exp 



rvr n 



(6.2) 



of [11] is fulfilled. We then choose a = ^ (that minimizes ac- 1 ) getting by |19J 
for k < e (Q " A) the equality (FTTH) and 

£ |$ 9 (x)|<Mr) (6.3) 

X:X(T)>1 

The invariance of the Q q under translations leads to (|3.14jl . On the other hand the 
last inequality gives 

E i^wi< E E ( 6 - 4 ) 

X:suppX9z T:supp T3x X:X(r)>l 

< E Mr)<^r<^- Q (6.5) 

T:supp T3x 

by choosing e~ A = ne~ a and taking into account that a — 1 > 1; here the first 
sum is over all multi-indexes X whose support contains a given point x. This 
implies (l3~T6l) and also J-TTfl since |<r(A | < |9A| J2x: mp p xbx \®q( X )\- Note 
furthermore that one easily gets the bound (|4.11j) from relations l|fi.lj) and (|6.3j) 

We shall also gives a bound needed below. Define the diameter of a contour F as 
diam T = max Xj3/g x(r) d(x, y) where d(x, y) is the distance between the site x and 
y. Then 

oo 

E I^WI< E ^r)<Y,^ Xn a- n <{ K e- a ) A (6.6) 

X:suppX3x r:S(r)93; 71>A 

diam supp X>A diam T>A 

We now turn to the proof of Theorem [2j We put h q = —p q and 

a q = h q — min h m (6.7) 

m 

The boundary condition g is called stable if a q = 0. Our first step is to show that 
if the boundary condition q is stable then all g-contours are stable implying that 
S^(A; fj,i, fi 2 ) coincides with S g (A; /i 2 ). 

We notice that when a q < 1, then the condition 

^ m (Int m r;/ii,/i 2 ) ^ e 2|aint m r| /g g N 

Hg(Int m r;/i 0)/ ui) ~~ 
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with m 7^ q for a g-contour T ensures that this contour is stable, provided 

Indeed by ((SUEj) flWx — 9q is bounded by a g + 2«;e Q . Since |<9lnt m T| < |L(r)| 
we get taking into account the Peierls estimate (13.71) : 

flmax|g(r)| -( flma x-fl g )|g(r)| -m(Int mixt £j /Xj, /i 2 ) 



z,(r) = w(r) e - Smax|iii)| e 



H g (Int m i x t T; /ii, /z 2 ) 

< g -(/3J-5)L(r) 

For a volume A, we use diamA = max^gA d(x, y) to denote its diameter The 
following proposition show that all g-contours are stable whenever a q = 0. 

Proposition 1. Assume (3 is large enough so that e~ /3J+5 = e~ a < ^d+i)*, ' ^ en 



(i) : if a m diamA < 1, and a q = 0, then 

^m(A; /ij, fl 2 ) < e a m |A|+2Ke- Q |9A| 

H 9 (A;/ii, fjv) 

(ii) : if a q = 0, then 

1 — m (A 

H g (A;/ii,/x 2 ) 

(iii) : if a m diamA < 1. then 



(6.9) 



< e 3«e-|9A| ( g 1Q) 



^ m (A;/ii,/x 2 ) < e (i+ 5Ke - Q )|9A| 
H m (A;/ii,// 2 ) 

'['lie proof is analog to that of Theorem 3.1 in |2Hj using our previous estimates. 
We give it below for the reader's convenience. 

We first introduce the notion of small and large contours. We say that a m- 
contour Y is small if a m diamT < 1; it is large if a m diamT > 1. We also define 
the partition function S™ all (A) which is obtained from H^(A) by replacing the sum 
over stable contours in (13. 9 J) by a sum over small contours. If we sum instead, only 
over contours which are at the same time small and stable, we denote the resulting 
partition function S' small (A). Finally we will use the shorthand notation S m (A) for 
1 — m (A 

We will show the three items of the proposition inductively on diamA. 

Thus we assume that (i), (ii), and (iii) have already been proved for all volumes 
with diam A < k. 

Proof of (i) for diamA = k 

For any contour V in A, and any in, we have diamlntmT < k — 1. We can use 
the inductive assumptions (ii) and (iii) that all g-contours and all m-contours are 
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stable. Therefore 

S ? (A) _ S' g (A) 



S m (A) ~' m (A) 



(6.12) 



Using the convergence of cluster expansion ()3.16|) and definition ()6.7|) , one imme- 
diately gets (i). 

Proof of (ii) for diam A = k 

To control the ratio S m (A)/S ? (A), we shall rewrite the partition function H m (A) 
using relation H3.1J3.2J) . Consider for a set of compatible m-contours in A, the family 
{Ti, . . . , r n }gxt ge of its mutually external large m-contours. The others contours are 
small m-contours in Ext = A \ \JiV(Ti) or any m-contour in Ext = UjExtTj. 
Therefore 

n 

S m (A) = Sr all (Ext)n^( r i) 5 m(Intr,) (6.13) 

r r r ilarge i=l 

Dividing and multiplying by [JlLi e 9max|5(ri)l S g (Int Tj) = S q (Int) n"=i e 9q ^ s ^ r ^ , we 
get 

S "*( A ) _ \p 5^ all (Ext)S g (Int) "rr g a \s(Tj)\ (r v -wStt^l S m(Int Tj) 
E q (A)~ ^ Eq (A) 11 ^(Intr,) 

y v ' rp p i large ^ v 7 «=1 ^ v ' 

U 1,— J 1 n/ext 

(6.14) 

Note that all g-contours in A and all small m-contours in A are stable by the 
inductive assumptions (i) and (iii) respectively. Therefore the various partition 
functions in the first factor of the right-hand side of (J6.14I) are equal to the cor- 
responding truncated partition functions, which can be controlled by convergent 
cluster expansion. We get by (J3.16I) 

5^ iaJ1 (Ext)5 g (Int) jj^jgto)! < e -ftr all |Ext|+^|A\int| JJ gjfeiscroi 

"S^) i=l i=l 

xe Ke- a (|9A| + |9Int| + |aExt|) 

where h s ™ aU is the free energy obtained from 5™ a11 . Using the facts that \h q +g q \ < 
Ke~ a , \S{T)\ < \L{T)\ and bounding |<9lnt| + |<9Ext| by |<9A| + 2d £\ L(r<), we find 



7 small 



(Ext)5 q (Int) -q ^,15(^)1 < e -(^ all -ft q )|Ext| e Ke-"[2|aA|+(l+2d) Ei LCTi)] (6.15) 



i=l 



Combining this bound with (I6.14JI . the Peierls estimates (13.71) . the inductive 
assumption (ii), and <9lnt T < 2dL(T), we get: 
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E m (A) 



(/c™»-V|Ext| 77 e -PJ\L(ri)\ e K8d+i)Ke-°>]L(ri) 



1—1 /r, r \ lar 8 e i=l 

n 

< e 2«e-«|c>A| e _ (/l small_ fe?) | Ext | JJ e _( a+ i) L ( ri ) ^g 16 ^ 

rr, r \ lar se i=l 

where for the last inequality we used the hypothesis (2c? + l)/te~ Q < 1 and a = 
f3J-5. 

At this point we need a technical lemma proved in [T2| (see the proof below) 
Lemma 1. Consider the partition function 

^ A )= E n^)e L(ri) (6.17) 

{Ti,...,r n } i=l 

of a gas of contours with activities 

Let —Is — — lim A | Z d (1/|A|) In Z(A) be the corresponding free energy. Then 



n 



II K ( T ^ ^ e Ke " a+1 l aA l (6.18) 

where the sum is over contours in A provided 

a>s (6.19) 

< ^ (6.20) 

To apply this lemma to ()6.16|) we put 

a = a + l (6.21) 

a = h 5 ™ 311 -h p = a m + hT n ~ h m (6.22) 

and 

^-{r ,r, !fr^3 ^ 

For Ke~ a+1 < 1, the Mayer expansion for In Z(A) is convergent. Using the fact 
that it only contains large contours (diamT > l/a m for each contour contributing 
to In 2(A)) one has (see (|6.6jl ): 

s < (Ke~ a )^ (6.24) 
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Moreover the difference ft-™ a11 — h m is the free energy of a gas of large contours with 
again diamT > l/a m and thus again for Ker a ~ x < 1 shares the same upper bound 

\h s ™ al1 - h m \ < (Ke~ a )^ (6.25) 

Hence 

a-s >a m -2{ne~ a )^ (6.26) 

Therefore the assumption a — s > will be fulfilled if a m > 2 (/te~ a ) am i.e. if 
(a m /2) am > ne~ a . This is actually true because the function (a/2) a has a minimum 
at a = 2/e for which it takes the value e~ 2//e ~ 0.47 that is greater than the upper 
bound l/(2d+ 1) required for ne~ a . 

Applying the lemma (with the value 1)6.21(1 ) to ()6.16)l immediately gives ()6.10j) . 

Proof of (in) for diam A = k 

The inequality ()6.11)) follows immediately from 1)6. 9|) . 1)6.10)) , and the fact that 

a m |A| < a m diamA \dA\ < \dA\ (6.27) 

Proof of LemmaU\For ne~ a+1 < 1, the partition function 2(A) can be controlled 
by convergent cluster expansion. In particular for the interior Int= Ujnt Ti of a set 
of external contours {Ti, T n } ext we have the estimate 

n 

Z(Int)e- ? l Int l > e - Ke " 5+1 l aint l > J] e -™«r B+1 \m)\ ( 6 . 28 ) 

i=i 

Here the first inequality stems from 1)3.17)1 and the last from the hypothesis 1)6.20)1 . 
Therefore 

n n 

e~ a \ Ext \ \{ z(T,) < e_a|EXt ' II ?(r,)e 2dKe " 5+lL ( r ^i(Int)e J l Int l 

{ri,-,r„} ext i=i {rx,...,^}^ i=i 



■i i 



< e -iAi n 5r ( r ^ e(2 * te " 5+l+s)L(ri) ^( int ) 

{r lv ..,r n } ext i=i 

< e- 5|A| Z(Int) 

where in the last inequality, we used s < Ke~ a+1 . This gives (16.18)1 using again the 
estimates (13.17)1 . □ 

This ends the proof of Proposition [TJ 

We now come back to the proof of Theorem[2j We put [i\(t) = fii+t, ^(t) = /i2+t 
with e w + e^ 2 = 1. Then, by definitions lETSj) . (l3~T71 . and (Q, we have 

a em p - ctmixt = JWt - Pemp = t + ln(e w + e^ 2 ) + lim -|- [lnZ mixt (A) - lnZ e (A)] 

A|Z d |A| 

(6.29) 
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where 

n 

Z q (A) = e-^ A l ~;(A) = E II 4( r *) ( 6 - 3 °) 

{ri,...,r n > j=i 

I j- j j us C omp 

The function t + ln(e Atl +e^ 2 ) = t is obviously increasing, negative for t < 0, positive 
for t > 0, and it intersects the horizontal coordinate axis only at one point t = 0. 
The difference a m ; xt — a emp will satisfy the same properties with the intersecting 
point slightly changed) provided 

1 



IM 



d d 
— In Zemp(A) - ^ In Z mixt (A) 



< 1 



(6.31) 



uniformly in A. 

Let us first give an upper bound on the derivative ^z' q (T) of the truncated ac- 
tivity. By virtue of relations (|3.8JI . (|3.6j) one gets for every stable g-contour T, that 
either -^z'JT) = 0, or: 

d , /-p\ d d H m (Int m T) 
_ln W (r)-^(TW_l 



dt 5 ? (Int m r) 

< (i + |int m r|)<(r)<|v(r)| e - a i L ( r )i 



<(r) 



On the other hand, 
1 d 



\A\dt 



lnZ ? (A) 



9_j 
dt 



4(r) 



X 



^(A\{r» 



^(A) 



(6.32) 
(6.33) 



<- £ 

Here the sums are over contours V containing a given point x and 

z q (A\{v}) = J2* 

{ri,...,r n } comp 

where the sum goes over all families {Ti, •••,r n } comp compatible with T. This sum 

can be bounded using the cluster expansion. Indeed, denoting S the set of sites at 
distance less or equal to 1 from the support of T, we get by (I3.13j) . (l3.17l) : 



(6.34) 



^(A\{r}) 

Z q {A) 



expJ- £ $ g (X)l<expj|Sl Y, 

< exp{(2d + l)\S(T)\Ke- a } < exp {(2d + l)L(l>e~ Q } (6.35) 

Inserting this bound in (|6.33j) and taking into account the inequality ()6.32|) and 
the estimate 



\V(T)\ < |5(r)|diamr < \S(T)\ 2 < e 



\S(T)\ 



< e 



L(T) 



gives 



1 <9 , t~r / A \ 

HS lnZ « (A) 



e -aL(r)+(2d+2)L(r)Xe- Q! 

T:supp Tga; 



< z/"e ~"e 

n>l 



n — an 2 (rf~ 1 ) kc a n 
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Using that K ~ 4.9z/, we get that for 2(d + l)ne < 1 this last sum is less than 
1/2 . This implies (16.31^ ending the proof of Theorem [21 ■ 

7. Proof of Theorem 3 
Let us first explain the idea of the proof. 

In dimension d = 3, the first excitation of a flat interface is obtained by replacing 
an empty site of the interface by an occupied one (denote 1^ such an interface) or 
by replacing an occupied site of the interface by an empty one (denote I^ own such 
an interface). In the first case, we have 

w (/) e - p W>f3)|S(i)| 

/:supp IC.V 
supp /=supp 1^ 

n=0 v 7 

„*„—5/3Ji I „*„-5/3J2 

= (c * e -M + C * 2 e-^) N C 4 +C f 



In the second case we have 

u(I)e-*M'WW\ 



J:supp ICV 
supp /-supp /down 



( w - 4 )e^i- /3Jl ) n e-^- /3j2 )( 7V - 4 - n ) /4 

n=0 m=l ^ 

( c * e -WJi +c ; e -2/3J 2 )4 



e 2m(^-/3Ji) e -(8-2m)(^-/3J 2 ) 



(c*e- pJ1 + c%e 



(c*e-^ + cte-P^Y 



In dimension d = 2, the first excitation of a flat interface is obtained by splitting 
it into a left part and a right part and then shift the right part by a height 1 or — 1 
(denote J e i such an interface). Then, we have 

J2 wWeTterimr)] 



/:supp ICV 
supp /—supp / e | 



^Z 1 / AT — m\ p {vl-l3-h)np-{V*2-PJ2){N-rn-l) 



ra=0 



* -2/3Ji I r * p -2f3J 2 

v 1 2 ; c\e~^ + cJje-^ 

Note, that in fact such modified interfaces are not allowed with the boundary 
condition ^ mixt > em P. However, one can lightly modify this boundary condition to 
allow such interface, leaving the resulting surface tension unchanged. 
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To study the difference F'(V) — Ffi &t (V), we shall express the quantity 

/:supp ICV 

as the partition function of a gas of excitations to be called walls or more generally 
aggregates with small activities at low enough temperature. 

This will allow us to exponentiate this quantity and to obtained that the differ- 
ence (F'(V)— Ffi at (V))/N can be expressed as a convergent series (up to a boundary 
term) whose leading terms are : 

1 c*e" 5/3Jl + c*e~ 5/3j2 1 (c*e^ Jl + c*e" 2/3j2 ) 4 
~P c\e-P J i + c* 2 e-^ ~ p (cje-^i + c* 2 e~^) A 
in three dimensions and 

2 c\e- 2f5Jl + c* 2 e~ 2l3J ' 2 
~P c\e-^ + c* 2 e~^ 

in two dimensions. 

7.1. Decorated interfaces. The first step is to express the quantity A(I) defined 
by (gj§ in a form suitable for our purpose. We write for each support of the mixt 
or emp cluster C in A(I): 

? frr . IcnDI j / s~f\ \CnD\ 

£ lr n \cnu\ s fr ,s \Cnu\ 

e~* emp(C) ^ = 1 + ( e ~ $emp(c) — - 1) = 1 +^ emp (C) 
Define for a connected family A of support of clusters: 

VWt(A) = Vw(C), V> emp (v4) = J"J 4 mp (A) 
ceA CeA 

Then 

cns(/)^0 cns(/)^0 

n m 
{Ai,...,A„} oomp :A i nS(/)^0 i=l {A 1 ,...,A m } comp :A i nS(/)^0 j=l 

where the sums are over compatible families of connected sets A{ of support of 
clusters touching the interface. As it was done for multi-indexes, it is convenient 
to sum all A with the same support say V to be called decoration. We define the 
weight 

n 

</VxtOP)= Vw(^4)= H^(d) (7.1) 

A:suppA=© {Ci,...,C n }:Ud=V i=l 

n 

4mp(^)= VW04)= Il^empCCi) (7.2) 

A:suppA=D {Ci,...,C„}:Ud=V i=l 
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This leads to 



-PF'{V) 



J2 w (j) e -*M*s>ra Ihw(A) 

/:SUPP ICV {'E'l, -.,Dn}comp: 1=1 

'D 4 nS(/)^0,£> i n.D#0 

m 

{"D 1 ,...,r> m } comp : J = l 



where the sums are over compatible families of decorations touching the interface. 

We define a decorated interface as a triplet J de = 2) m ixt) ®em P }, where / is an 
interface, 2? m i x t is a collection of mixt-decorations touching the interface and S) emp 
is a collection of emp-decorations touching the interface. 

_§ . cc)l CnD 

The weights of decorations may be controlled with the inequalities \e mixU ; i c i 

-1| < (e-l)|$ mixt (C)| and \ e ~^ v) V -1\ < (e - l)|$ emp (X?)|. Together with 
(|4.1ip . this implies the bounds (see [29]): 

|^ mix t(^)| < (8e(e - l)Ke- a ) L ^ (7.4) 
IVW(2>)| < ( 8e (e - l)ne~ a ) L ^ (7.5) 



7.2. Walls and aggregates. We now introduce the notion of walls and aggregates 
by the following definitions. 

Consider a decorated interface J de = {/, 2D mixt , 2) emp }. Let Il denote the hori- 
zontal hyper-plane Xd = and let 7r denote the projection parallel to the vertical 
axis on this hyper-plane: 7r(xi, xa) = (xi, Xd-i)> 7r (^4) = ^xeA^(x). A pair 
{$x,Sy} of the interface I is called correct if 

(1) tt(x) = ir(y) (the bond xy is vertical). 

(2) there is no other pair pairs {s x >, s y >} of the interface such that ir(x') = n(x). 

(3) there are no decorations such that it(T>) D ir(x). 

The connected components of the set of correct pairs are called ceilings and denoted 
C. The connected components of the set of non correct pairs are called walls (see 
Fig. 5). 

For a wall W, we use S a (W) to denote the set of sites for which s x = a. The 
set supp W = S Q (W) U Si (WO U S 2 (W) is called support of the wall W. As for the 
contours and the interfaces, we will also use S(W) = Si (W) U 52 (W) to denote the 
set of occupied sites of the wall, L\(W) (respectively L 2 (W)) to denote the number 
of nearest neighbour pairs (x,y) such that such that s x — 1 and s y = (respectively 
s x = 2 and s y = 0) and L{W) = L X {W) + L 2 {W). 

Denote further W(J) the set of walls of the interface /. 
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The union Uiyew(/)supp W U£, g/ de V split into maximal connected components 
called aggregates. Namely, an aggregate w of a decorated interface is a family 

w = {W u ...,W n ;V u ...,V m } 

such that the set U" =1 supp Wi UTLi Vj is connected. 

A collection of aggregates {w±, . . . , w n } is called admissible if there exists a dec- 
orated interface J de such that Wi, . . . , w n are the aggregates of 7 de . 

Here comes a difference for definitions between dimension 2 and dimension 3. 

In 3-dimensions an aggregate is called a standard aggregate, (or aggregate in 
standard position), if there exits a decorated interface J de such that w is the unique 
aggregate of J de . To any aggregate w, we can associate a unique standard aggregate 
which is just a translate of w. In this way, one can associate to any admissible 
collection of aggregates a unique collection of standard aggregates. Such collections 
are called admissible collections of standard aggregates. To an admissible collection 
of standard aggregates we can associate in a unique way an admissible collection 
of aggregates. 




2 11 








2 1 




2 
2 
2 
1 
1 









1 

2 
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Figure 6: The walls corresponding to the interface of Figure 4. 

7.3. Expansions. With these definitions and notations, one gets from ()7.3|) . the 
following expansion 



{w 1 ,...,w n } adln i=l 



(7.6) 



where the sum is over admissible collections of standard aggregates and the activi- 
ties of aggregates are given by 

e -pj 1 Li(W)-pj 2 L 2 (W)+ f j,* 1 \S 1 (W)\+ t i%\S 2 (W)\ 

z(w) = [[ eP{ ^ 2) \ s{w) \^ e _ pJl - c * e -f3j 2 yn(w)\ 11 ( 7 - 7 ) 



vew 



Here ip{T>) = ^Vxt^) for the mixt-decorations and ip{T>) = V'emp(^) for the emp- 
decorations. 
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In two dimensions we proceed differently. For an aggregate w, we consider the 
ceiling components to the left and to the right of w respectively denoted and 
C^ R \ Let (respectively h^) be the second coordinate of any empty site of 
(respectively C^). We define the position p(w) of th aggregate w as p(w) = 
and the height of the aggregate w as h{w) = — h^ L \ An aggregate w is now 
called standard aggregate (or aggregate in standard position) if p(w) = 0. One 
obviously can associate to any aggregate an aggregate in standard position and to 
any admissible collection of aggregates a unique collection of standard aggregates. 
As before such collections are called admissible collections of standard aggregates 
and to an admissible collection of standard aggregates we can associate a unique 
admissible collection of aggregates. 

For an admissible collection of aggregates {w'±, . . . , u^jadm, one has the constraint 
Y^i=i h( w 'i) = 0- If we l e t, for i = 1, . . . , n, Wi denote the standard aggregates 
corresponding to w'^ this constraint reads also Y^i=i h( w i) = 0- From (|7.8j) . we 
then get: 

n n 
{u;i,...,-u)„} adm i=l i=l 

where the activities z{w) of aggregates are also given by (|7.7|) . 

We now introduce the notion of elementary walls which are the walls correspond- 
ing to the elementary interfaces mentioned in the beginning of this section. 

In two dimensions, an elementary wall W e \ is a wall that contains two pairs 
{s x ,s y } and {s x ,s z } such that the site x is occupied, the sites y and z are empty, 
y is above x and z is to the left or to the right of x. 

The activity of such a wall are easily computed and we have 

? z ( w ei) = \ ri ; «t (7-9) 

W el :suppW el ={ X ,y,z} 

In three dimensions, we call elementary wall, either a wall (corresponding to the 
elementary interface J"j P ) that contains four pairs {s x , s Xi }, i = 1, . . . ., 4, such that 
the site x is occupied, the sites Xi are empty and all the sites leave on the plane 
0, or a wall (corresponding to the elementary interface /down) that contains 



3 



X 

eight pairs {s x , s Xi }, {s Xi , s yi } , i = 1, , 4, such that the sites x and Z{ are empty, 

the sites yi are occupied, the sites x and Xi leaves on the plane x 3 = — 1, the sites 
yi leaves on the plane x 3 = 0. 

In the first case we have 

? Z M = * Rl * ~R1 — 7.10 

W el :suppW el ={x, X i, X 2,X3, X 4} 
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while in the second case 

( c *e- 2 ^ Jl + c *e- 2/3j2 ) 4 



E 



z{w e i 



W el :suppW el ={x,x 1 ,x 2 ,x 3 ,X4,yi ,2/2,2/3,2/4} 



(7.11) 



For non elementary walls and aggregates, the activities 1)7.7)1 may be bounded as 



z(w) < 



(a g e-P J ) L ^ 



V z(w) < — 



{a g e-P J ) L 



(7.12) 
(7.13) 



(In 



w:L{w)=L V 1 2 7 

where 

Be(e - l)Ke -a e /w if > 2d 
1 otherwise 

This allows to exponentiate the partitions functions of the gas of aggregates in 
the RHS of (|7.6|) and 1)7.8)1 . for low enough temperatures. We define multi-indexes 
Y corresponding to aggregates as function from the set of aggregates into the set 
of non negative integers, and we let supp Y = U^y^^isupp w. The truncated 
functional corresponding to the activities z is given by 

a(Y) ■n*H y( " ) ( ? - 14 ) 



U w Y(w)\ 

where the factor a(Y) is defined as in 1)3.12)1 . Then, 



E Ylz(wi)=exp 

{wi,...,w n } aAm 1=1 



{ E 



-/3iV^+<T(V|*) 



(7.15) 



Y:supp ycv 



in 3-dimensions and 

n 



E n^KE^ ) = pr{E^) = } ex p{ E 

{wi,.. .,«)„} adm j=i i=i i=i y :S upp ycv 

n 

= Pr I E] ^(^j. 



Y:supp YCV 

Q \ e -fiNF+*<ym (716) 



in 2-dimensions. Here 



E 

Y:suppY90 



*00 



E 

Y:supp Ynv c ^9 



supp K n n 1 

supp Y n ir(y)\ 



and 



Pr{$>( 







|supp Y n n | 
E { ^ 1 ,...^ n } adffl EKLi ^K)<H EILi MO, 



(7.17) 
(7.18) 



i=l 



{l«l,...,W)„} adl 



nr=i z ( w i) 



(7.19) 
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The series (|7.17l) converges whenever 8e(e — l)K 2 e~ a < 1. This can be seen by 
verifying the convergence condition (|6.2j) with the activities z(w) and with the 
contours replaced by aggregates. 

Notice that the multi-indexes involved in the sum of (I7.18|) intersects both ^{V) 
and V = 7L d \ V. Thus at low temperatures, this sum can be bounded by L d ~ 2 
times a constant so that this term will give no contribution to the surface tension 
in the thermodynamic limit. In addition, the probability Pr{^™ =1 h(iUi) = 0}, may 
be controlled by known techniques (see e.g. [HI HSU EH E2|) and 

n 

lim (1/N) lnPr{V/i(wi) = 0} = 

i=l 

Hence, by taking the thermodynamic limit L — > oo in equations ()7.6j) and ()7.8|) we 
obtain: 

e -/3(w, emp -^) = c * e -PJi _|_ c * e -P J 2 (7.20) 

Whenever a multi-index Y contains only one aggregate w (Y(w) = 1 and Y(w') = 
for w' 7^ w) one has *f?(Y) = z{w). The leading terms of T are then obtained with 
the help of relations H7.9jl - H7.lljl . This gives the expressions Q5.5J) and Q5.6jl . Taking 
furthermore into account the expression H4.16j) ends the proof of the theorem. 
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leculaire, Universite de Mons-Hainaut for warm hospitality and financial support. 

References 

[1] E. A. Guggenheim. Trans. Faraday Soc, 41:150, 1945. 

[2] R. Defay and I, Prigogine. Trans. Faraday Soc, 46:199, 1950. 

[3] J. G. Eberhart. J. Phys. Chem., 70:1183, 1966. 

[4] B. von Szyszkowsky. J. Phys. Chem., 64:385, 1908. 

[5] H. P. Meissner and A. S. Michaels. Ind. Eng. Chem., 42:2782, 1949. 

[6] A. W. Adamson. Physical Chemistry of Surfaces. John Wiley and Sons, New- York, 1997. 
[7] R. Defay, I. Prigogine, A. Bellmans, and D. H Everett. Surface tensions and absorption. 

Longmans, Green and Co., London, 1966. 
[8] J. De Coninck and J. Ruiz. Interfacial tensions for binary mixtures versus Solid-On-Solid 

models. Preprint, 2003. 
[9] J. C. Wheeler and B. Widom. J. Chem. Phys., 52:5334, 1970. 
[10] J. L. Lebowitz and G. Gallavotti. Phases transitions in binary lattice gases. J. Math. Phys., 

12:1129, 1971. 

[11] G. Gallavotti. Phase separation line in the two-dimensional Ising model. Commun. Math. 
Phys., 27:103, 1972. 

[12] R. L. Dobrushin. Gibbs states describing the coexistence of phases for a three-dimensional 

Ising model. Theory Prob. Appl., 17:582, 1972. 
[13] M. Blume, V. Emery, and R. B. Griffiths. Phys. Rev. A, 4:1071, 1971. 

[14] D. Ruelle. Statistical Mechanics: Rigorous Results. Benjamin, New- York, Amsterdam, 1969. 
[15] D. Ruelle. Thermodynamic Formalism. Addison Wesley, Reading, 1978. 
[16] Ya. G. Sinai. Theory of Phase Transitions: Rigorous Results. Pergamon Press, London, 1982. 
[17] M. Zahradnik. An alternate version of Pirogov-Sinai theory. Commun. Math. Phys., 93:359, 
1984. 



ON THE STATISTICAL MECHANICS AND SURFACE TENSIONS OF BINARY MIXTURES 37 

[18] A. Martin Lof, G. Gallavotti, and S. Miracle-Sole. Some problem connected with the co- 
existence of phase in the Ising model. In Statistical mechanics and mathematical problems, 
Lecture Notes in Physics, volume 20, page 162. Springer, Berlin, 1973. 

[19] S. Miracle-Sole. On the convergence of cluster expansion. Physica, A 279:244, 2000. 

[20] R. Kotecky, and D. Preiss. Cluster expansion for abstract polymer models. Commun. Math. 
Phys., 103:491, 1986. 

[21] J. Bricmont, J. L. Lebowitz, and C.-E. Pfister. On the surface tensions of lattice systems. 

Ann. N. Y. Acad. Sci., 337:214, 1980. 
[22] J. Bricmont, J. L. Lebowitz, and C.-E. Pfister. Non-translation invariant Gibbs states with 

coexisting phases III. Commun. Math. Phys., 66:267, 1979. 
[23] S. Miracle-Sole. Surface tension, step free energy, and facets in the equilibrium crystal. J. 

Stat. Phys., 79:183, 1995. 
[24] D. B. Abraham. Surface structures and phase transitions-Exact results. In Phase Transitions 

and Critical Phenomena, volume 10. C. Domb and J. L. Lebowitz, editors, Academic Press, 

London, 1986. 

[25] C.-E. Pfister. Interface and surface tension in Ising model. In Scaling and Self- similarity in 

Physics, J. Frohlich, editor. Birkauser, 1982. 
[26] H. N. V. Temperley. Proc. roy. Soc. A, 199:3, 1949. 
[27] Work in progress. 

[28] C. Borgs and J. Imbrie. A unified approach to phase diagrams in fields theory and statistical 

mechanics. Commun. Math. Phys., 123:305, 1989. 
[29] C. Dobrovolny, L. Laanait, and J. Ruiz. Surface transitions of the semi-infinite Potts model 

I: the high bulk temperature regime. J. Stat. Phys., 114:1269, 2004. 
[30] J. Bricmont and J. Frohlich. Statistical mechanical methods in particle structure analysis of 

lattice field theories II. Commun. Math. Phys., 98:553, 1985. 
[31] R. L. Dobrushin, R. Kotecky, and S. B. Shlosman. Wulff Construction: A Global Shape from 

Local Interactions. Providence, New- York, 1992. 
[32] J. De Coninck, R. Kotecky, L. Laanait, and J. Ruiz. SOS approximants for Potts crystal 

shapes. Physica, A 189:616, 1992. 



